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The  paper  discusses  the  potential  incorporation  of  higher  performance  gyroscope 
and  accelerometers  in  the  horizontal  channels  of  Litton's  local-level  inertial 
positioning  system  and  the  resulting  improvements  in  positioning  and  deflection 
of  the  vertical  determination.  It  further  provides  deflection  error  estimates 
based  on  the  use  of  simple,  and  advanced  data  reduction  methods  and  indentifies 
limitations  of  interest.  '  1 
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Ci  POTENTIAL  HARDWARE  AND  SOFTWARE  IMPROVEMENTS  OF 

ZD  INERTIAL  POSITIONING  AND  GRAVITY  VECTOR  DETERMINATION* 
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rH 

< 


H,  BAUSSUS  VON  LUETZOW 


ZUSAMMENFASSUNG 

In  dlesem  Vortrag  werden  der  potentlelle  Elnschluss  von  Krelseln  und 
Beschleunlgungsmessern  hoeherer  Genauigkelt  in  den  horlzontalen  Kanaelen 
des  vertikal  orientlerten  Inertialortsbestimmungssystems  von  Litton 
und  resultierende  Verbesserungen  hlnslchtllch  Posltlons-und  Lotabwelch- 
ungsennittlung  eroertert.  Ferner  werden  Lotabweichungs-Fehlerschaetz- 
ungen  mittels  einfacher  und  hoeherer  Datenreduzierungsverfahren 
geliefert  und  interessferende  Beschraenkungen  Identifiziert. 


SUMMARY 

The  paper  discusses  the  potential  incorporation  of  higher  performance 
gyroscopes  and  accelerometers  in  the  horizontal  channels  of  Litton's 
local-level  Inertial  positioning  system  and  the  resulting  improvements 
in  positioning  and  deflection  of  the  vertical  determination.  It 
further  provides  deflection  error  estimates  based  on  the  use  of  simple 
and  advanced  data  reduction  methods  and  Identifies  limitations  of 
interest. 

RESUME 


Cette  etude  est  une  discussion  sur  1 'incorporation  potentlelle  des 
gyroscopes  de  haut  rendement  et  des  accelerometres  dans  les  voles 
horizontalles  du  systeme  Inertlelle  de  determination  du  point  a  niveau 
locale  de  la  compagnle  Litton.  Cette  communication  fournit  une 
evaluation  de  l'erreur  de  deflection  basees  sur  l'emploi  de  m§thodes 
fondamentalles  et  avancSes  de  reduction  de  donnees.  Les  limites  de 
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TRODUCTION 

e  present  Rapid  Geodetic  Survey  System  (RGSS),  developed  by  Litton 
stems,  Inc.  for  the  U.S.  Army  Engineer  Topographic  Laboratories  (ETL), 
s  the  capability  of  determining  horizontal  and  vertical  positions, 
flections  of  the  vertical,  and  gravity  anomalies  with  average  rms 
rors  of  1  m,  0.3m,  Earcsec,  and  2  mgal,  respectively,  for  50  km 
ns  under  utilization  of  post-mission  adjustments.  It  operates  as 
quasi  local-level  system.  It  thus  does. not  require  altitude 
mplng,  permits  Kalman  stochastic  enrol*  control  without  great 
mplexlty  under  consideration  of  observed  velocity  errors  at 

—  - *Tnv+t 

ed  paper,  to  be  presented  on  Aug.  17,  1981. 
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vehicle  stops,  .and  employs  effective  post-mission  adjustments 
with  the  aid  of  terminal  position  and  gravity  vector  Information 
only.  Present  critical  hardware  consists  of  an  A-1000  vertical 
accelerometer,  two  A-200  horizontal  accelerometers,  and  two 
6-300  gyroscopes.  Although  data  utilization  from  repeated  runs 
has  resulted  In  average  deflection  rms  errors  of  about  1  arcsec, 
the  requirement  of  a  maximum  rms  error  of  0.5  arcsec  necessitates 
the  Incorporation  of  higher  performance  gyroscopes  and  acceler¬ 
ometers.  This  would  simultaneously  achieve  conventional  surveying 
accuracy  (1 0~5 ) . 

2.  SIGNIFICANCE  OF  HIGHER  PERFORMANCE  GYROSCOPES  AND 
ACCELEROMETERS  FOR  POSITION  AND  VERTICAL  DEFLECTION 
ERROR  REDUCTION 

Following  the  approximate  elimination  of  the  effects  of  constant 
gyro  biases  on  positions  under  utilization  of  accurate  Initial 
and  terminal  coordinates,  the  remaining  dominant  error  source 
is  gyro  correlated  random  noise.  In  this  respect,  the  pertinent 
parameters  for  the  6-300  instruments  are  a  standard  deviation 
of  0.002°  hr-1  and  a  correlation  time  of  3  hours.  Another  error 
source  to  be  considered  is  the  accelerometer  scale  factor, 
assumed  to  be  constant  for  a  test  run.  The  standard  scale 
factor  error  for  the  A-200  accelerometer  Is  0.01%.  Correlated 
accelerometer  noise  Is  characterized  by  a  standard  deviation  of 
10  mgal  and  a  correlation  time  of  40  minutes.  Both  correlated  gyro 
random  drift  and  correlated  accelerometer  measurement  errors 
affect  the  accuracy  of  deflections  of  the  vertical  estimated 
from  Inertial  data  and  Initial  and  terminal  deflection  components, 
possibly  augmented  by  corresponding  geodetic  azimuths.  In 
order  to  achieve  deflection  accuracies  smaller  than  or  equal 
to  0.5  arcsec  rms.  It  is,  therefore,  necessary  to  Install  higher 
performance  gyros  and  accelerometers.  A  significant  reduction 
of  the  two  autocorrelation  parameters  is  also  a  necessary 
prerequisite  for  the  optimal  statistical  estimation  of  deflections 
of  the  vertical.  The  Litton  6-1200  gyro  Is  expected  to  be 
compatible  with  the  stringent  deflection  accuracy  requirement. 

Its  random  drift  error  is  about  O.OOl^r-1  and  its  autocorrelation 
parameter  Is  short.  For  details,  reference  Is  made  to 
Litton  [1974]  and  Huddle  [1977a].The  A-1000  accelerometer 
addressed  by  Litton  [1973,1975]  has  a  standard  scale  factor  error 
of  approximately  0.005%.  Correlated  accelerometer  noise  Is 
2  mgai  rms,  equivalent  to  aboutO.35  arcsec.  As  mentioned  In 
the  Introduction,  this  Instrument  Is  already  Incorporated  In  the 
vertical  channel  of  the  present  RGSS.  In  conjunction  herewith, 
an  Improved  velocity  quantizer  has  been  Installed.  These 
Improvements  have  resulted  in  altitude  errors  of  0.3m  rms 
and  gravity  errors  of  2  mgal  rms  after  2- hour  runs  and  under 
utilization  of  post-mission  adjustments.  Of  further 
significance, the  precision  of  the  6-300  gyros  and  A-200 
accelerometers,  primary  horizontal  channel  components  of  the 
present  R6SS,  Is  greater  than  Indicated  by  the  accuracy 
parameters.  This  has  been  observed  by  Schwarz  [1979*1  and 
by  ETL  after  evaluation  of  repeated  test  runs  during  1980. 
Multiple  runs  can  also  be  expected  to  reduce  the  single  run 
standard  deviations  associated  with  an  advanced  RGSS  (AR6SS), 
advocated  by  the  author  in  1978.  * 
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3.  BASIC  METHOD  FOR  THE  DETERMINATION  OF  DEFLECTIONS 
OF  THE  VERTICAL 

A  presentation  of  both  Kalman  error  control  and  post-mission 
adjustment  and  of  the  determination  of  deflection  components  has 
been  made  by  Baussus  von  Luetzow  [1979],  With  respect  to  position 
determination,  Schwarz  [1979b]  noted  that  the  errors  Increased 
considerably  during  tests  In  strongly  mountainous  terrain.  This 
Is  due  to  the  fact  that  under  these  conditions  gravity  anomalies 
and  deflections  of  the  vertical  cannot  be  satisfactorily  modeled 
as  quasl-statlonary  random  variables  In  the  real  time  Kalman 
error  controller.  Fortunately,  these  Inherent  difficulties  do 
not  significantly  affect  the  determination  of  deflections  of 
the  vertical, the  detailed  discussion  of  which  requires  the 
utilization  of  the  local-level  error  differential  equations, 
applicable  to  horizontal  motion  and  restricted  to  land  vehicles: 


d  - 

siy  * 

-*y 

(1) 

i-x  , 

lx 

dt 

P. 

(2) 

iL  y  = 

dt  9 

SN*Z  -  +6^  +  aE 

(3) 

d  • 

—  x  « 
dt 

"  sEez  +  g$E  -  gt  +  aN 

(A) 

• 

y  2  * 

tn  i-  (fijj  +  pjjsec  *)  x  +  u 

(5) 

ar*»- 

• 

R  +  wz*E+  B 

(6) 

Jr*" 

• 

+  wz»k  +  Vz+  y 

(7) 

For  simplicity,  the  symbol  6  In  front  of  the  dependent  variables 
has  been  ommltted.  The  applicable  coordinate  system  Is  evident 
from  figure  1 . 

x  (north) 


A 


>  y  (east) 


Figure  1 

Applicable  Coordinate  System 


W  ■ 
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Symbols  used  In  the  foregoing  equations  are  with  z  downward: 

x  north  angular  position  error 

y  east  angular  position  error  (for  longitudes  west  of 
Greenwich) 

x  north  velocity  error 

y  east  velocity  error 

a  azimuth  axis  angular  drift  rate  error 

B  north  axis  angular  drift  rate  error 

y  east  axis  angular  drift  rate  error 

g  normal  gravity 

f  geographic  latitude 

R  mean  earth  radius 

®z  azimuth  platform  attitude  error 

»N  platform  tilt  error  about  north  axis 

*£  platform  tilt  error  about  east  axis 

gn  product  of  g  and  deflection  component  n 

g£  product  of  g  and  deflection  component  & 

dVv 

V— a  north  acceleration  of  survey  vehicle 

dVF 

S  « — £  east  acceleration  of  survey  vehicle 
E  dt 

aE  correlated  east  accelerometer  error 

aN  correlated  north  accelerometer  error 

nN  -  ncos4>  north  earth  rate 

pn.Ve/R  north  angular  rate 

“N*nN  +  pm  ■  Hcos4  +  Ve/R  north  spatial  rate 

«£ ■  pE  *  -V„/R  east  spatial  rate 

+  p2  »  flsln*  4  ve/R  •  tn$  vertical  spatial  rate 

In  Inertial  land  navigation  terms  involving  »E  in  equations  (5) 

and  (6)  are  neglected  In  Lltton's  Kalman  error  controller.  The 
Initial  conditions  at  time  to  ■  0  are 
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■  0,  * 


(0) 

N 


*>0* 


>(0).  x(0).  y(0).  0 


In  the  general  RGSS  node  of  operation,  position  and  gravity  vector 
Information  Is  generated.  The  corresponding  Kalman  mechanization 
Is  then  simplified  by  setting  ■  u»  *  0  In  equations  (5)-(7) 
which  results  In  the  uncoupling  of  tne  horizontal  channel  error 
differential  equations.  This  Is  not  critical  for  travel  periods 
of  2  hours. 

As  an  example,  equations  (4)  and  (7)  may  then  be  formulated  at 
the  end  of  the  first  travel  period,  whence  t  ■  ti. 


* 


-  *i  ♦  aNi 

-  XI  +  f  vdt  ♦  e0 
r  J 


(8) 

(9) 


In  eq.  (8),  the  existence  of  an  Initial  accelerometer  measure¬ 
ment  error  has  been  taken  Into  consideration. 

The  magnitude  Xi  Is  well  estimated  by  the  observed  velocity  error 
xi.  After  Implementation  of  the  tilt  correction  (est.),  equations 
(9)  and  (8)  read 

•i}'1  “  -«i  +  J  \ dt  +  C0  (10) 


*1°  «  g(-«i  -  g^i  -  C0\  +  aNl  -  aN()  (11) 


The  estimation  error  Is  generally  negligible. 

Because  of  the  quasl-ellmlnatlon  of  Integrated  terms  Involving 
A,  y%  x  In  equations  (5)-(7)  by  appropriate  tilt  corrections, 
constant  gyro  biases  may  be  obtained  with  a  high  degree  of  „ 
approximation  by  a  solution  of  these  equations  with  x  »  y  «  x  *  0 
and  o  ■  5,  6  ■  I*  y  •  Y»  It  Is  then,  according  to  Huddle  [1977bJ, 

*E  "  *(*  +  ®  Y*)  (12) 

£*  t  +  B  ^  (13) 

'■■'('T*  *  ‘  ) 


(14) 
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Under  consideration  of  terminal  deflection  and  azimuth  closure 
errors.  It  Is  possible  to  compute  a,  8 ,  7  from 


*  *  estimated  "  *  observed 


n  estimated  *  n  observed 


^  estimated  ^  observed  (^) 


where  A  designates  geodetic  azimuth. 

Examination  of  equations  (12)-(14)  reveals  that  In  moderate 
latitudes  there  exists  relatively  small  coupling  for  time  Intervals 
not  exceeding  2  hours.  For  this  reason,  linear  approximations 
♦e«*  yt  and  8t  have  been  used  with  success.  Highly  accurate 
deflection  determinations,  accomplished  by  means  of  an  ARGSS,  would, 
however,  require  consideration  of  a  terminal  azimuth  error. 

For  t  *  tj,  equation  (11)  may  be  formulated  with  *  0  as 

x(ti)  «  x±  -  g?El  -  9  Ui  -  i0)  +  g  /  1  (v-7)dt  +  aNi  -  (18) 


i(tn)  -  xn  «  9®^  -  9  Un  -  6e)  ♦  9  /  (y-T)dt  +  #Nn  -  aNo  (19) 

o 


Equations  (12)  and  (19)  show  that  the  computation  of  *Ei  is  to 
the  first  order  associated  with  an  error 

5?Ei  *  -  r  C  9  /tn(Y‘7)dt  +  aNn  *  aN03  (20) 

XI  o 


The  basic  solution  for  the  prime  deflection  of  the  vertical  from 
equations  (18)  and  (19)  Is 


(21) 


Equation  (21)  may  be  supplemented  by 


«o  +  £ 


(6C-  -  -  (*.  *  71  *e  ) 

e  °  si  tn  Sn 


(22) 


to  account  for  astrogeodetlc  deflection  errors  and  accelerometer 
bias  errors.  For  a  straight  traverse,  the  last  term  In  eq.  (22) 
tends  to  cancel  out. 

The  rms  deflection  error  orfc^)  can  be  computed  by  covariance 
analysis  Involving  the  terms  without  parentheses  in  eq.  (21). 

Under  Inclusion  of  the  first  two  terms  of  eq.  (22)  It  Is 

var  *  var  4  var  y*  +  0  -  |^)2  var  £0  4  (51)2  var  (23) 


where  var  ai  Is  the  accelerometer- Induced  variance  and  var 
designates  the  gyro-induced  variance.  Under  consideration  of  present 
RGSS  parameters.  Identified  in  section  2,  the  rms  deflection  error 
c£.  is  approximately  presented  In  Figure  2. 


Distance  Traveled  (km) 


Figure  2 

Approximate  RMS  Deflection  Error  as  Function  of  Travel  Time 


The  variation  of  normal  gravity  g  with  altitude  Is  not  critical 
In  terrestrial  applications.  The  scale  factor-generated  error 
does  not  vanish  with  respect  to  L-shaped  traverses.  It  Is, 
however,  relatively  small  for  travel  times  not  exceeding  2  hours. 
Heading  sensitivity  Induced  by  significant  azimuth  changes 
appears  to  be  somewhat  critical  and  may  require  empirical 
corrections.  Schwarz  [1979b]  discussed  the  Impact  of  heading 
sensitivity  on  positioning  and  concluded  that  error  reductions  may 
be  achieved  by  a  modified  smoothing  procedure.  The  problem  Is  more 
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Intricate  In  the  case  of  deflection  determinations.  The  use  of 
static  accelerometer  measurements  In  conjunction  with  Initial 
and  terminal  deflections  together  with  a  simplified  Kalman 
filter  originated  with  Huddle  [1977a3.  The  above  analytical 
presentation  Is  one  of  several  quasl-llnearizatlon  techniques 
according  to  the  terminology  of  Schwarz  [19803.  It  may  be 
refined  by  a  finite  difference  solution  of  the  system 
differential  equations  (l)-(7)  with  +  •  const,  with  or  without 
periodic  filter  corrections.  This  would  result  In  weight  factor 
expressions  as  a  function  of  time  involving  the  generating 
variables  e  •  C[x(t),  y(t)3,  n  •  n[x(t),  y(t)3,  »E  '  «E  •  ~  *n  • 

a,  6,  and  y.  Approximate  constant  vehicle  velocities  would  f 
simplify  the  computations.  A  matrix  solution  involving  all 
unknown  c  s  and  n’s  would  then  replace  the  independent  single 
component  solutions. 

4.  ADVANCED  METHOD  FOR  DEFLECTION  DETERMINATION 

Equation  (21)  under  Inclusion  of  the  first  two  terms  of  eq.  (22) 
or  a  comparable  refined  solution  Identified  In  section  3  may  be 
formulated  as 

Ej  *  it  *  \  (24) 

mm 

where  Ei  *  +  T£l  -  ^  Is  a  message  variable, 

Is  a  signal  variable,  and  the  remaining  terms  denote  a  noise  variable  -  nlt 
The  collocation  method  in  physical  geodesy  permits  In  semi-flat 
terrain  the  estimation 

Ee  *  E  «i(Ei  +  nt)  «  Ajfij  4  Nt)  (25) 

where  Is  the  matrix  of  regression  coefficients  at  to  be  computed. 

It  Is  then  in  matrix  form,  with  bars  Indicating  covariances, 

+  *k>  -  0,  1,  ...n  (26) 


The  solution  for  the  regression  coefficient  matrix  follows  as 

♦  W’1  <27> 


Due  to  the  sizable  gyro  and  accelerometer  correlation  times  associated 
with  the  present  RGS5,  the  advanced  method  characterized  by  equations 
(25)  and  (27)  does  not  provide  significantly  better  deflection 
estimates.  Considerably  Improved  estimates  would,  however,  result 
under  utilization  of  Intermediate  deflection  constraints.  In  the 
context  of  an  ARGSS,  l.e.,  under  consideration  of  very  short  correla¬ 
tion  times,  the  advanced  method  may  be  expected  to  yield  better 
estimates.  In  addition.  It  constitutes  the  statistical  framework  for 
an  optimal  area  solution  under  utilization  of  data  relating  to 
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several  traverses.  The  computation  of  topographic  corrections  and 
of  spatial  covariance  functions  for  applications  In  strongly 
mountainous  terrain  has  been  presented  by  Baussus  von  Luetzow  [1981], 
Alternatively,  It  would  be  necessary  to  conduct  repeated  runs.  In 
this  way,  rms  deflection  errors  not  to  exceed  0.35  arcsec  are  within 
the  reach  of  an  ARGSS. 

5.  CONCLUSION 

The  Incorporation  of  Identified  or  other  higher  performance  gyros  and 
accelerometers  and  of  Improved  velocity  quantizers  in  Litton *s  RGSS 
would  result  in  an  ARGSS.  Non-linear  gyro  bias  corrections, 
elimination  of  errors  caused  by  gyro  heading  sensitivity,  and 
advanced  data  reduction  techniques  or  repeated  runs  are  then  expected 
to  facilitate  deflection  determinations  with  standard  errors  not 
to  exceed  0.35  arcsec  for  a  travel  time  of  2  hours  or  about  60  km 
distance.  As  a  by-product,  conventional  surveying  accuracy  would 
be  achieved.  Deflection  changes  £-£0  and  n-n0  with  standard  errors 
of  the  order  0.1  arcsec  for  the  same  travel  time  and  vehicle  speed 
would  require  supplementation  by  a  high-performance  gravity 
gradiometer. 
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